UCRHEP-T506 
RECAPP-HRI-2011-003 

Signatures of supersymmetry with non-universal Higgs mass at 

the Large Hadron Collider 



Subhaditya Bhattacharya^" 1 , Sanjoy Biswas* 2 , Biswarup Mukhopadhyaya* 3 
<N ■ and Mihoko M. Nojiri* 4 



^Department of Physics and Astronomy 
University of California, Riverside, California 92521, USA 

^Regional Centre for Accelerator-based Particle Physics 
Harish- Chandra Research Institute 
SpH 1 Chhatnag Road, Jhunsi, Allahabad -211 019, India 

c-| ■ * Institute for the Physics and Mathematics of the Universe 

University of Tokyo, Chiba 277-8583, Japan 
^ ', Theory Group, KEK, 1-1 Oho, Tsukuba, Ibaraki 305-0801, Japan 

■ The Graduate University for Advanced Studies (SOKENDAI) 

^ \ 1-1 Oho, Tsukuba, Ibaraki 305-0801, Japan 

in 
o 

Abstract 

We discuss large non-universality in the Higgs sector at high scale in super symmetric 
theories, in the context of the Large Hadron Collider (LHC). In particular, we note 
that if rriH u 2 — m H d 2 is large and negative (~ 10 6 GeV 2 ) at high scale, the lighter 

C$ • 

slepton mass eigenstates at the electroweak scale are mostly left chiral, in contrast 
to a minimal supergravity (mSUGRA) scenario. We use this feature to distinguish 
between non-universal Higgs masses (NUHM) and mSUGRA by two methods. First, 
we study final states with same-sign ditaus. We find that an asymmetry parameter 
reflecting the polarization of the taus provides a notable distinction. In addition, we 
study a charge asymmetry in the jet-lepton invariant mass distribution, arising from 
decay chains of left-chiral squarks leading to leptons of the first two families, which 
sets apart an NUHM scenario of the above kind. 
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1 Introduction 



With the Large Hadron Collider (LHC) already running, one feels closer than ever to glimpses 
of physics beyond the standard model (SM). Supersymmetry (SUSY) [1, 2] has always re- 
mained an attractive hunting ground in this context. The LHC has brought added impetus 
to not only the search for SUSY, but also the more ambitious proposal to identify the over- 
seeing high-scale physics that can lead to typical low-energy spectra. Such high-scale physics 
is often envisioned as the 'organizing principle' behind the plethora of low-scale parameters 
in the minimal supersymmetric standard model (MSSM) [3, 4] to be seen at low energy. A 
frequently adopted approach in this direction is to embed MSSM in a minimal supergravity 
(mSUGRA) scenario [5] , where all the low-scale parameters can be generated from: 

m i/2, m o, A), tan/3 and sgn(ii) 

where mi/2,mo, Aq are the universal gaugino mass, scalar mass and trilinear scalar coupling 
parameters respectively at the high-scale, tan/3 = (H u )/(Hd) is the ratio of the two Higgs 
vacuum expectation values and /i is the SUSY-conserving Higgsino mass parameter in the 
MSSM superpotential. 

However, the mSUGRA model can be branded over-simplistic, as the assumption of 
universality doesn't follow from any known symmetry principle. For example, gaugino mass 
non-universality can occur in supersymmetric Grand Unified Theories (SUSY-GUT) [6] with 
non-trivial gauge kinetic functions [7, 8]. Non-universality in the scalar sector can also be 
motivated from the 5*0(10) D-terms [9], apart from the phenomenological requirement to 
keep CP-violation and flavour changing neutral currents (FCNC) under control [10]. Also, 
in mSUGRA, one assumes that the Higgs mass parameters have their origin in the same 
m Q which generate squark and slepton masses, which is completely ad hoc. For example, in 
SUSY-GUT theories based on the 50(10) group, sfermions and Higgs fields belong to differ- 
ent representations and can therefore arise from independent high-scale mass parameters. In 
view of this, one can, within the SUGRA scenario itself, expect the Higgs mass parameters to 
arise from high scale value(s) different from m . Thus models with non-universal Higgs mass 
(NUHM) are of considerable interest, and their viability in respect of both collider signals 
and issues such as the dark matter content of the universe has been recently investigated 
[11, 12, 13]. 

One can incorporate the non-universality in the Higgs sector in two different ways. In the 
first kind, one can have both of the soft Higgs mass parameters originating in a high-scale 
value m' which is different from mo, the universal high-scale mass for squarks and sleptons. 
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In other words, one can postulate rn 2 Hu = rn 2 Hd = m' ^ m\ [12]. On the other hand, it is 
also possible to have H u and Hd evolve down from two different high-scale inputs. In the 
later case, the high-scale SUSY parameters are given by: 

m 1/2 ,m ,m 2 Hu ,m 2 Hd ,A ,t&n[3 and sgn(fx) 

The split between the two Higgs squared masses at high scale introduces additional fea- 
tures in the running of various mass parameters down to the electroweak scale. In its most 
drastic manifestation, such a situation can give rise to the sneutrino (z>) as the lightest super- 
partner of standard model particles. Since a sneutrino dark matter candidate is disfavoured 
from available results on direct search, one then has to postulate the sneutrino(s) to be 
the next-to-lightest supersymmetric particle(s) (NLSP), and, for example, gravitino as the 
lightest supersymmetric particle (LSP). With this achieved, most of the allowed region of 
the NUHM parameter space leads to the right amount of relic density [14]. 

In this work, we propose using the LHC data to distinguish those cases where the super- 
particle spectrum in NUHM is most strikingly different from the usual mSUGRA scenario. 
As we shall see in the next section, this happens for a large negative high-scale value of 
mil 2 — m H d 2 - It not only leads to a large splitting between the left and right chiral sleptons, 
but also leads to the lighter slepton mass eigenstate of any flavour being dominated by the 
left chiral component. 

This feature, marking a drastic departure from the expectations in mSUGRA, can be 
reflected in the signals of staus through the polarization of the taus that are produced either 
in their decay or in association with them [15, 16, 17]. In addition, the above hierarchy 
between left-and right-chiral sleptons can be probed by studying the spin correlation of jets 
and leptons produced in cascade decays of squarks. This correlation, as we shall see, affects 
the angular distribution of the lepton in x° — * I l T , manifested through certain measurable 
kinematical variables [18, 19, 20]. 

To explain further, the large splitting between the left-and right-chiral sleptons sometimes 
yields a hierarchy where the right-chiral ones become much heavier than not only the left- 
chiral ones but also the low- lying chargino/ second lightest neutralinos over a large region 
of the NUHM parameter space. Thus they are hardly produced in collider experiments. 
At the same time, the (dominantly) left-chiral stau and the corresponding sneutrino being 
considerably lighter — even lighter than the lightest neutralino — the taus produced in their 
association are dominantly left-handed. This is due to the fact that the gauge couplings 
involved in the decay are chirality conserving, so long as one has large gaugino components 
in the lighter neutralinos and charginos. 
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The consequences that we focus on are two-fold. First, one notices the practically ubiqui- 
tous r in SUSY signals. Secondly, the signals often bear the stamp of left-polarized r"'s, in 
the products of their one-prong decay. With this in view, we analyze the polarization of the 
taus produced in the SUSY cascades in the same-sign di-tau (SSDr) final states associated 
with hard jets and missing transverse energy (fir)- We show how this leads to noticeable 
differences between the NUHM and mSUGRA spectra in the LHC environment. 

Furthermore, we study the polarization dependence of the angular distribution of the 
lepton produced in \% decay, which shows up in the charge asymmetry in the m q i distribution. 
Though the effect tends to wash out due to the presence of antisquark decay, nevertheless it 
can be observed at the LHC as more squarks are produced than antisquarks. 

This paper is organized as follows. We discuss various aspects of the model under con- 
sideration in the following section and identify the region of the m — m^j parameter space 
where the lighter stau is dominantly left-chiral. As we shall see below, this is achieved for 
large negative values of S. We choose a few benchmark points for our collider simulation. 
Tau-polarization and its implications are discussed in section 3, while the analysis revealing 
the chirality information on sleptons of the first two families is outlined in section 4. The 
numerical results for each of the two analyses mentioned above, based on a simulation for 
the 14 TeV run of the LHC, is presented in section 5. We summarise and conclude in section 
6. 



2 Features of the NUHM scenario and our choice of 
benchmark points 

2.1 Salient features of the scenario 

We consider the general case of NUHM, having a two-parameter extension of the mSUGRA 
scenario, in which the soft SUSY breaking masses m 2 Hu and m 2 Hd are inputs at high scale. 
The most important thing to remember here is that the renormalisation group evolution 
(RGE) of soft scalar masses is in general modified by the presence of a non-zero boundary 
value of the quantity S, defined as [3] 

S = m 2 Hu - m 2 Hd + Tr [m 2 Q - m 2 L - 2m* + + m|] (1) 
We assume universality in the sfermion masses, so that S = — m 2 Hd is high scale 
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boundary condition. The running of soft scalar masses of the third family squarks and 
sleptons are given at the one-loop level by [3] 

dm 2 Q3 2 



l^M 2 - 3g 2 2 M 2 2 - fgjM 2 + ±g 2 S + y 2 X t + y 2 X b ) (2) 

(3) 

y/M - fgjM! - 2 -g\S + 2y 2 X t ) , (4) 

dm- 2/4 1fi 1 \ 

' -glMl-^glMl+ V -glS + 2ytxX (5) 



dt 16tt 



d m i R 2 / 16 16 2 



dt 16tt 2 V 15 ai 1 3 " 6 6 5' 

^ = A, f-\ 2 M 2 - 3 5 2 M 2 - *-&S + yiX T ). ((3) 
dt 16tt 2 V 5 yi 1 y2 2 10 yi 

dm| fl 2/12 



A - 5^ + ^ + 2^). (7) 

where the notations for squark, slepton and gaugino masses have their usual meaning, and 
t = log(<5), yt,b,r are the t, b and r Yukawa couplings, and 

X t = ml 3+ m? R+ ml u+ A 2 t , (8) 
X b = ml 3 +ml R + ml d + Al (9) 
X r = m 2 L3 +ml R + m 2 Hd + A 2 T (10) 

Mass parameters of the first two family scalars run in a similar manner, excepting that 
the Yukawa contributions are vanishingly small. The main difference in the SUSY particle 
spectrum with respect to an mSUGRA scenario is the non-vanishing boundary value of S. 
If this boundary value is large in magnitude, the effect on the spectrum at low scale is 
naturally a rather pronounced departure from mSUGRA. Since the contribution of the term 
containing S comes with different factors in the running of left-handed squarks (sleptons) 
and right-handed squarks (sleptons), due to different U(l) hypercharge assignments, one can 
have large splitting in the left-right sector within each generation when j^l is substantially 
large. 

One can see from equation (6) and (7) that the effect of non-universal Higgs mass is 
rather pronounced in the slepton sector, the primary reason being that the running masses 
are not controlled by the strong sector. The most important difference it makes to the 
spectrum is that, for large negative values of S ( O(TeV) 2 ) ) [11, 12], the left-chiral sleptons 
tend to become considerably lighter than their right-chiral counterparts. This is in striking 
contrast to both mSUGRA and gauge mediated SUSY breaking (GMSB). An immediate 
temptation that the phenomenologist faces, therefore, is to extract some signature of this 
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'chirality swap' in the lightest sleptons at the LHC, which may put a distinctive stamp of 
NUHM on them. This, of course, has to be done with the help of leptons that are produced 
either in association with the low-lying sleptons or in their decays. Since the helicity of 
leptons of the first two families is difficult to measure in the collider environment, we feel 
that it is our best bet to latch on to the copious number of taus arising from SUSY cascades, 
and concentrate on those features of their decay products that tell us about their helicities. 

As has been noted already, the above effect is seen for large negative S. Such values of S 
therefore become the benchmarks for testing the special features of NUHM, and it is likely 
that in such condition only its footprints are noticeable at the LHC. Thus we examine next 
the kinds of spectra ensuing from large negative S, and look for their observable signature. 

A large negative S at high scale affects the running of the third family SU(2) doublet 
slepton (both the stau and the tau-sneutrino) masses in the same way as is done by their 
Yukawa couplings, thus bringing them down substantially at low energy. As a consequence, 
one can have both of them of the same order as, or lighter than, the lightest neutralino (x?)- 
In the latter situation, the left-chiral tau-sneutrino is lighter than the corresponding stau 
due to the SU(2) breaking D-terms (for tan/3 > 1) : 



m\ L = m|-cos(2/3)m|(--sin 2 M (11) 
ml = m| + cos(2/3)m|.i (12) 

In such cases, the tau-sneutrino has to be the NLSP, due to its unsuitability as a dark 
matter candidate as laid down by direct search results. A gravitino, for example, can be 
envisioned as the LSP and dark matter candidate in such cases. The lighter stau mass 
eigenstate can in principle also become the NLSP through mixing of the left and right chiral 
fields. However, this happens only in very restricted regions of the parameter space, as large 
mixing requires tan f3 to be on the higher side, a feature that is highly restricted in NUHM by 
the requirements of absence of tachyonic states as well as of electroweak symmetry breaking. 

The role of S in the running of m 2 Hu and m% is described by 

3$M + ^gfS + 3/*X t ) , (13) 
3g 2 2 M 2 - ^gfS + 3f 2 X b + / r 2 X r ) , (14) 

One can see above that a negative S tends to partially cancel the effects of top quark 
Yukawa coupling in the running of m 2 H and make it positive at low energy. m 2 H , on the 
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other hand, is routinely rendered positive at low scale due to the gauge interactions, and 
the effects of the term proportional to S often fails to make it negative as one comes down 
to the electroweak scale. Consequently, radiative electroweak symmetry breaking at the 
right energy requires a negative value of m 2 Hu at high scale. Of course, one is led to have a 
sufficiently large magnitude of \i to ensure that m? Hu + jj? remains positive at high energy. 

2.2 The choice of benchmark points 

As has been already explained, our purpose is to suggest some observations at the LHC, 
which will bring out the distinctive characteristics of the NUHM spectrum. Such distinction 
is most pronounced when the chiralities of the low-lying sleptons are reversed with respect 
to the corresponding cases in mSUGRA. This, we have found, is best achieved (and one 
is indeed optimistic about clear distinction) when S is large and negative (~ 10 6 GeV 2 ). 
For smaller magnitudes of S (< 10 5 GeV 2 ), the f L component of fx decreases, and the 
collider signature of this scenario is relatively less distinct. With this in view, the region 
in the parameter space with more than 90% of ft in fx has been shown in Figure 1. This 
region offers the best hope for recognising NUHM if SUSY is detected at the LHC. We 
have accordingly chosen some benchmark points for the study reported in the subsequent 
sections. Out of the regions answering to our chosen criterion, we have selected points with 
three possible mass hierarchies: 

rfiv Th < tti x o < m fl 
m 9TL < m fl < m x o 
m x o < m i>TL < m fl 

Our benchmark points (BP) NUHM-1 - NUHM-3 (shown in Table 1) are taken from 
three regions of the parameter space, corresponding to each of the above hierarchies. The 
code SuSpect (version 2.41) [21] has been used for this purpose. Two-loop renormalisation 
group equations have been used for running the mass parameters down to low energy, with 
the default option (namely, \J t^) for the electroweak symmetry breaking scale. The spectra 
are consistent with low energy constraints [22, 23] such as those coming from b — > S7 and 
the muon anomalous magnetic moment, and also with those from LEP-2 limits, such as 
m x ± > 103.5 GeV, mj± > 98.8 GeV and m h > 111 GeV. Electroweak symmetry breaking in 
a consistent fashion has been taken as a necessary condition in the allowed parameter space. 
For the case with Xi LSP, the requirement of relic density consistent with the recent data 
has also been taken into account in choosing the benchmark point (s) [24]. 
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Benchmark points 


NUHM-1 


NUHM-2 


NUHM-3 


Input 


m = 300 


m = 80 


m = 300 


parameters 


m i/2 = 300 


m i/2 = 460 


m 1/2 = 280 




tan/3 = 10 


tan p = 10 


tan/3 = 7 




170 


154 


154 




552 


437 


551 




151 


132 


133 


' L 


119 


106 


116 


m fl 


139 


124 


137 


m f2 


537 


424 


543 


m v o 

Al 


120 


187 


112 


m v o 


234 


361 


216 


mo 

A3 


939 


982 


950 




944 


987 


954 




234 


361 


217 


m xt 


944 


987 


955 


m~g 


734 


1066 


691 




645 


826 


618 


m h 


814 


1018 


791 


m d L 


741 


986 


706 


m d R 


742 


962 


710 




737 


984 


701 


m u R 


591 


879 


549 




111 


114 


112 



Table 1: Proposed benchmark points for the study of the NUHM scenario with m 2 Hu = —1.10 x 
10 6 GeV 2 and rnfi d = 2.78 x 10 6 GeV 2 . All the mass parameters are given in units of GeV. The 
value of Aq is taken to be zero and sign of \i to be positive for all of the benchmark points. 
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Benchmark points 


C 1 T T /~1 T"> A 1 


«~ CTT r~* T"> A o 


mbUGKA-3 


Input 


m = 80 


m = 350 


m = 300 


parameters 


mi/2 = 250 


mi/2 = 300 


mi/2 = 350 




tan/3 = 40 


tan ft = 40 


tan/3 = 10 


m e~ L i m fiL 


389 


362 


284 


m e R , m jl R 


363 


322 


202 




381 


354 


271 


mo TL 


353 


329 


269 


m fl 


283 


238 


197 


m f2 


377 


358 


285 


m v o 


99 


120 


140 


m v o 


183 


224 


261 




333 


392 


455 


m v o 


352 


409 


473 




182 


224 


262 


m x 2 ± 


353 


410 


473 


m 5 


b2S 


726 


8ol 


^t! 


AQA 


525 


573 


m h 


615 


683 


764 


m d L 


654 


720 


776 




636 


697 


745 




649 


716 


771 


m u R 


636 


698 


747 




111 


112 


111 



Table 2: mSUGRA benchmark points obtained based on similar cross-section in the same-sign 
ditau channel (mSUGRA-1 and mSUGRA-2) and in the opposite-sign same-flavor dilepton channel 
(mSUGRA-3). All the mass parameters are given in units of GeV. The value of Aq is taken to be 
zero and sign of fj, to be positive for all of the benchmark points. 
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We have obtained the mSUGRA BP's for comparison with the NUHM points using the 
criterion based on similar event rates (within ±30% tolerance) in two different channels. For 
the case where the distinction between these two scenarios is done using tau-polarisation, 
we have compared the event rates in the same-sign ditau (SSDt) channel for choosing 
our mSUGRA points. For the analysis based upon lepton-charge asymmetry, the event 
rates in the opposite-sign same-flavor dilepton (OSSFDi) channel have been compared as 
a benchmarking criterion. All the three NUHM BP's have been used for the first case and 
for the second case, only NUHM-1 and NUHM-3 have been considered, as the hierarchies 
mentioned above are not relevant for analysis based on lepton-charge asymmetry. Thus, we 
have obtained mSUGRA-1 which corresponds to both NUHM-1 and NUHM-3 following the 
criterion mentioned above in the SSDt channel and mSUGRA-3 corresponds to NUHM-3 in 
the OSSFDi channel. The benchmark point mSUGRA-2 corresponds to NUHM-2 having 
similar rates in the SSDt channel. Two of the chosen mSUGRA points (mSUGRA-1 and 
mSUGRA-2) are approximately compatible with the observed relic density. 

The values of various SUSY parameters in the chosen points are listed in Table 1. One has 
to further assume in the case of v T -NLSP and gravitino (G) LSP that the decay u T — > v T G 
does not a lifetime exceeding the age of the universe. The gravitino mass has to have 
accordingly allowed values, as dictated by the hidden sector of the overseeing theory [25]. 

3 Tau polarisation 

The signal of left-polarised tau is expected to be a very good discriminator between scenarios 
with NUHM and its universal counter part. Tau lepton plays a crucial role in the search for 
new physics. In particular, information on the chirality of a tau can be extracted following 
some standard procedures. The fact that the tau decays within the detector, in contrast 
to the electron or the muon, enables us to know about its chirality from the kinematic 
distribution of the decay products. In the massless limit where the tau is boosted in the 
laboratory frame, tau decay products are nearly collinear with the parent tau. In this limit, 
hadronic tau decays produce narrow jets of low multiplicity, to be identified as tau-jets. 
From the angle of polarisation studies, it is most cost-effective to work with the one-prong 
hadronic decay modes of the tau, which comprise 80% of its hadronic decay width and about 
50% of its total decay width. The main channels here are: 

T~ — > 7T~Z/ T 
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300 



200 400 600 800 1000 1200 1400 
m in GeV 

Figure 1: The allowed region for NUHM in the tuq — m\/2 plane for tanj3 = 10, m 2 Hu = —1.0 x 
10 6 GeV 2 , m \ i = 2.0 x 10 6 GeV 2 and Ao = 0. The light blue region is disallowed due to tachyonic 
stau and/or non-compliance of electroweak symmetry breaking conditions. The region on the left 
of the dashed line is also disallowed by constraints from b — > sj. In the region marked by +, one 
has m o < ra^ < rrif , whereas the pink region corresponds to m„ T < mo < mf The dark blue 
region has the hierarchy m„ TL < mf x < m^o . The lighter stau has 90% or more of tl over the 
entire allowed region. 

T~ — > (p~V r ) — )■ Tt~ll G V T 
T~ — > (di^r) — > 7r~7T 7r°Z/ T 

where we shall often denote both the p~ and the a{ by v . 

The first step in the extraction of polarisation information is to express some differential 

decay distributions of the r~ in the laboratory frame. Let the polarisation information 

be denoted by P T , where P T = ±1 correspond to taus with positive and negative helicity. 

Next, it is worthwhile to examine the laboratory frame variable z, defined as z = E ntV /E T , 

the fraction of the tau energy carried by the product meson. This variable can be related 

to 6, the angle between the direction of motion of the outgoing it~ or v~ and the axis of 

polarisation of the tau, which is taken to be along the direction of the tau momentum in the 

laboratory frame. In the limit E T » m T , 

„ 2z - 1 - c 2 

cosfl = — — 15 

1 — & 

where c = m v /m T . The expression for the case where the tau decays to the pion and a v T is 
obtained by setting m v = above. 
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The decay distributions in z for a r in the laboratory frame are given by [26] 

1 dT„ 



^ = [1 + P T (2z - 1)] (16) 



1 cftV m 2 m 2 r m 2 . 2 2 

' ! -[— §sm 2 u; + 1 + cos 2 w + P T cosO 



where 



T v dz (m 2 — m 2 ) (m 2 + 2m 2 ) m 2 

2 

( — £ sin 2 oj sin 2u tan# — 1 — cos 2 a;)], (17) 

ml m v 

1 dr„ T _ m 2 m 2 cos 2 w + sin 2 a; + P cos# 

r„ (m 2 — m 2 )(m 2 + 2m 2 ) m 2 

2 

772 777 

( — ?rCos 2 wH sin 2a; tan # — sin 2 a;)] (18) 

m 2 v m v 

(m\ - ml) + {m 2 T + m 2 v ) cos 6 

COSUJ = r^~> — 2 n , / 2 ^ — a ( 19 > 

(m; + m%) + (m; - m*J cos 6> 
In the experiment, one looks for hard jets from the tau, which corresponds to large values 
of z. A close inspection of Equations (16), (17), (18) shows that the energy distribution of 
the decay products from the decay of (P T = — 1) are in significant contrast to that from 
t r (Pt — +1) • When P T = +1, the hard r-jet consist largely of either a single pion or 
longitudinally polarised vector mesons iy£)- For P T = —1, on the contrary, the hard r-jet 
mostly comprises transversely polarised vector mesons only (i>t)- This conclusion becomes 
almost self-evident in, for example, the extreme case of collinearity, with cosw = 1 and 
sin a; = 0. 

It should, however, be remembered that the quantity z is not amenable to actual mea- 
surement in the detector, and therefore the distinctions pointed out above are still somewhat 
theoretical in nature. It is therefore necessary to translate the distinction in terms of mea- 
surable quantities. The energy distribution among the pions arising from the decay of the 
p~ and a>i offer such a variable. It is the variable R = E n /E p , the fraction of the energy 
of v carried by the charged pion. For the case where the p~ is produced in r~ decay, the 
distribution in R in the laboratory frame is given by [26] 

^1 . m 



dT(p L -» 2tt) 
dR 



(2P-1) 2 (21) 
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The distribution for a{~ is more complicated but has similar qualitative features. The 
reader is referred to [27] for the detailed expressions. The broad indication is that transversely 
polarised vector mesons favour even sharing of its momentum among the decay pions whereas 
longitudinally polarised ones favour uneven sharing of momentum among its decay products. 
Since the polarisation of the parent tau governs the level of polarisation of either type in the 
vector mesons v, the distribution in the variable R therefore is a reflection of the helicity of 
the tau whose signal one is concerned with. 

Obviously, one always has R = 1 when the tau decays as t~ — > tt~i> t . What one must 
utilise, therefore, is the difference in R-distributions between the cases with vt and v^. When 
the decaying tau has p T = +1, one should mostly have vl in the hard jets, in addition to 
the inconsequential single pions, giving its characteristic distribution on R. A contrast can 
be seen in the decay of a tau with p T = -1, where the hard tau-jets can be expected to be 
largely vt, with a different distribution in R. 

Hence, one can use the charged-pion spectra arising from the two-stage decays 



to probe the polarisation of the parent tau. We utilise this possibility to identify the NUHM 
spectrum in cases the low-lying stau is of left chirality, which attaches similar chirality 
(same as helicity at high energy) to the taus either arising from stau-decay or produced 
in association with it. With this in view, we have selected tau-jets in our simulation with 
Pt > 40GeV and \rj\ < 2.5, assuming a tau jet identification efficiency of 50%, with a fake 
tau jet rejection factor of 100 [16]. 

4 Lepton charge asymmetry 

Another discriminator which is sensitive to the mass hierarchy between the right-and left- 
chiral sleptons is the charge asymmetry in the jet-lepton invariant mass distribution [18, 19, 
20]. In the NUHM scenario (with large negative S), the lighter slepton mass eigenstate is 
dominated by the left-chiral component (Zi ~ II). Hence, for < m x o (a criterion mostly 
satisfied by the 'extreme' NUHM scenario considered by us), the leptons produced in the 
decay of x° will be mostly left-handed. In the usual mSUGRA scenario, on the other hand, 
one expects the leptons to be mostly right-handed as the lighter slepton mass eigenstate is 
dominantly right-chiral (li ~ Ir) and the decay proceeds via the Bino component of x\- 
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This feature can be exploited to unmask NUHM by studying the charge asymmetry in 
the lepton-jet invariant mass (m^) distribution produced in the squark decay chains, where 
li stands for the lepton reproduced in X2 decay. We shall consider sleptons of the first two 
generations only, for which left-right mixing is negligible, and the coupling of the leptons to 
the Higgsino components of a neutralino is also very small. 

In this section we describe the spin correlation in the following decay chain 



where I2 denotes the lepton produced in the subsequent step of the cascade. Due to the 
chiral structure of the squark-quark-neutralino coupling, the quark produced in the squark 
decay will be left-handed in the massless limit. The x° produced in decay is also polarized 
having the same helicity as that of the quark as they are produced from the decay of a scalar. 

In the rest frame of the squark produced in the initial hard scattering, a negatively 
charged lepton produced in the subsequent decay of the x° wm appear back-to-back or in the 
same direction as that of the quark depending on whether the slepton is left-chiral or right- 
chiral. 1 Exactly the opposite directional preferences hold for a (positively charged) antilepton 
vis-a-vis the quark produced in the chain. Therefore, we expect an asymmetry between 
the distributions m j7 - and Tn^+. This can be utilised to define the following asymmetry 
parameter: 



where % stands for the ith bin. A measurement Ai should thus yield information on the 
chirality of the low-lying slepton produced in the chain. 

However, there are some experimental difficulties involved in the measurement of such 
an asymmetry- 

1. In the decay of a the asymmetry in the lepton-jet invariant mass distribution has 
a sign opposite to that of the corresponding q L . This is because the left antisquark 
decays via gaugino coupling into a right-handed antiquark. Since jets initiated by a 
quark or an anti-quark are indistinguishable, it is impossible to disentangle the squark 
and antisquark production channels. However, the LHC is a pp machine where more 

1 The other inputs that go into this argument are (a) The x° produced in squark decay is sufficiently 
boosted, and (b) the X2 decays largely in the s-wave. 



(22) 




(23) 
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squarks are produced than anti-sqaurks, a significant 'net' charge asymmetry in the 
rrijh distribution can finally survive. All one needs in order to measure this charge 
asymmetry is a substantial excess in the production of q^g and q^q^ over pairs 
containing squarks and antisquarks, and also gluino pairs. 

2. In an experiment, it is not always possible to distinguish between the lepton (li) out 
of a X2 an d the lepton (Z 2 ) coming from slepton decay. We have taken the invariant 
mass distribution using the harder of the two leptons, a role in which li fits in most of 
the time. 

In NUHM, one expects negative charge asymmetries, whereas in the usual mSUGRA 
scenario they are expected to be positive, especially in the high invariant mass bins. However, 
in mSUGRA, depending on the mass hierarchy, the leptons produced in x° decay can also 
be dominantly left-handed if < < m x o, as the diagonal component ((6^)22) of the 
neutralino mixing matrix wins over (6^)21- m that case, one would expect a dip in the 
asymmetry distribution at a lower value of rriji and a peaking behaviour at the higher end. 
This is expected because the splitting between m x o and II is smaller than that between m x o 
and Ir. One can use this feature to separate an mSUGRA-type scenario. 

5 Collider simulation and numerical results 

We have simulated events for ^/s = 14 TeV, including initial-and final-state radiation, mul- 
tiple scattering etc. We have used parton distribution functions CTEQL6L1 [28] for our 
analysis, with the renormalisation and factorisation scales set at the average mass of the 
final state particles. 

5.1 Simulation strategy: ditau final states 

To study the polarisation of the tau in SUSY cascade for both NUHM and mSUGRA sce- 
nario we have used the code TAUOLA (version 2.9) [29] interfaced with the event generator 
PYTHIA (version 6.4.16) [30]. The spectrum has been generated using SuSpect (version 
2.41) [21]. TAUOLA has been suitably modified to incorporate the probability of producing 
left-or right-handed tau in the decay of SUSY particles. For cases where the v T and/or 
the f is lighter than the lightest neutralino, decay branching fractions of the lightest neu- 
tralino have been calculated using SDECAY (version 1.3b) [31] and fed into Pythia. The 
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finite detector resolutions have been taken into account following the specifications listed, 
for example, in [32]. 

The final state that we have considered is a pair of same-sign ditaus (SSDr), together 
with at least three hard central jets and large missing E T . Same-sign ditaus are preferred 
because they are less beset with SM backgrounds. We consider events where the taus have 
one-prong hadronic decays. 

The following cuts have been imposed on each event- 

• p T > 40 GeV, \rj\ < 2.5 for each tau jet. 

• pt > 100, 100, 50 GeV, \r)\ < 2.5 for the three associated jets, in decreasing order of 
hardness. 

• E T > 150 GeV. 

It should be reiterated that our main purpose is to obtain the observable difference 
between the NUHM scenario under consideration and an mSUGRA scenario. Situations in 
mSUGRA leading to tau-rich final states are most likely to fake NUHM phenomenology. 
Therefore, we have followed the criteria already mentioned in section 2.2, and isolated the 
regions where the total rate of SSDt+ > 3 jets + fir is within ±30% of the rate predicted 
for corresponding NUHM benchmark point. 

5.2 Simulation strategy: lepton charge asymmetry 

The charge asymmetry in the lepton-jet invariant mass distribution has been studied using 
the event generator HERWIG (version 6.5) [33] which takes into account the spin correlation 
in SUSY cascades. Spectra have been generated using IS A JET (version 7.78) [34] and the 
input parameters have been tuned in such a way that the spectrum generated is similar 
to that produced by SuSpect. A fast detector simulation has been done using AcerDET 
(version 1.0) [35] for reconstructing the isolated leptons, jets and fix, which also takes into 
account the finite detector resolution of the visible momenta. 

The final state under consideration is consists of a pair of isolated leptons of opposite 
charge and same flavor (OSSF) with more than three jets and missing Et, i.e., e + e _ + 
> 3 jets + fix- 

The preselection cuts [36, 37] imposed in this case are the following- 

• pr h > 20 GeV and p Th > 10 GeV, \rj\ < 2.5 for the two leptons. 
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• pt > 100, 50, 50 GeV, \rj\ < 2.5 for the three associated jets, in decreasing order of 
hardness. 

• M eff > 600 GeV where, M eff = fi T + S|p T | 

where, the summation is taken over all visible particles. 

• £ T >0.2M e// 

The SUSY backgrounds come mainly from two independent xt decay. One can eliminate 
this by taking the flavor subtracted combination e + e~ + ji + ji~ — e ± ji T and this cancels out 
the background contribution from the charginos up to statistical fluctuations. The Standard 
Model background, already small after imposing the above cuts, undergo further suppression 
in this process [37]. 

The leptons are combined with each of the two hardest jets and, for identifying the 
desired decay chain, the combination for which the invariant mass is smaller has been 

selected . The rriji± distribution for this subsample, for both the hard and soft lepton have 
been calculated. Depending on the mass splitting between the neutralino and slepton one of 
these leptons will be dominated by the 'correct' lepton, i.e., the one adjacent to the quark in 
the decay chain and will give the desired charge asymmetry in the jet-lepton invariant mass 
distribution. 

5.3 Numerical results 

Ditau final states: 

We first present the numerical results of our analysis using of the polarisation properties of 
the tau. In Table 2, we have tabulated the event rates for all the NUHM and the potentially 
faking mSUGRA points for the S'S'Dr-channel. Event rates have been predicted for an 
integrated luminosity of 100 fb~ l . After applying all the cuts to suppress the SM background, 
one has similar event rates for both the NUHM and corresponding mSUGRA points, which 
is not surprising because we have identified the mSUGRA points following the criterion of 
similar event rate. 

For the benchmark point NUHM-1, v Th is the LSP, and the lighter f\ is dominantly left- 
chiral. Taus are mainly produced in the decay of x\ ~^ rf(20.5%), xt ~* tz> tl (26.5%) and 
t — > r%5(100%) and therefore the taus are mostly left-handed. The contributions from x%i 
x\i xf an d ^2 are negligible as they are heavier in the spectrum. For NUHM-2 we similarly 
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NUHM-1 


NUHM-2 


NUHM-3 


mSUGRA-1 


mSUGRA-2 


31 


51 


28 


41 


46 



Table 3: Number of events in the SSDt channel at an integrated luminosity of 100 fb~ l after 
applying the cuts listed in Section 5.1, in addition to a cut on the R variable (R > 0.2 ^ for all of 
our benchmark points . 

have lighter stau mass eigenstate dominated by the left-chiral component but here the mass 
hierarchy between the f\ and the Xi is opposite to that of NUHM-1, i.e. m.f < m x o. At this 
benchmark point, Xi decays into 11 pair as well as vv pair including the third generation. 
The mass difference between the lighter stau and tau-sneutrino is less than m w , hence the 
decay proceed mainly via the two body decay mode ff* 1 — > i>*7r ± and the three body decay 
fi —¥ v^V^v^. However, final states with higher pion multiplicities also have non-zero 
branching fractions, but we have not taken into account these modes, as they do not change 
our conclusion. In NUHM-3, the LSP is the lightest neutralino, however we still have a light 
enough v TL . The lighter f l5 of course, dominantly left-chiral here. The taus produced in 
SUSY cascade therefore are mostly left-chiral for all the NUHM points. The corresponding 
R-distributions (taking into account the SM contributions) for the respective benchmark 
points have been shown in Figure 2. Thus the distinction criterion set down by us is seen to 
survive the washouts caused by various extraneous SUSY cascades. 





NUHM-1 


NUHM-2 


NUHM-3 


mSUGRA-1 


mSUGRA-2 


< 0.8) 


0.77 


0.74 


0.76 


0.71 


0.72 


2 (R > 0.8) 


0.23 


0.22 


0.24 


0.29 


0.28 


_ 0i(.R<O.8) 
O 2 (R<0.8) 


3.35 


3.36 


3.17 


2.45 


2.57 



Table 4: The ratio r for NUHM and corresponding mSUGRA scenario. 

In the corresponding mSUGRA benchmark points (mSUGRA-1 and mSUGRA-2), the 
lighter stau is dominantly right-chiral. However both the stau are heavier than the second 
lightest neutralino and lightest chargino. Taus are produced mainly via the decay of W 
and Z produced in the decay of X ° 4 (xt W*), (xlh/Z), X f -+ (x^W±), ( X f h/Z) and 
xt ~^ x\W ± . Hence the contributions to SSDt channel come from two same sign W-decay 
produced in SUSY cascade or one from W-decay and one in Z decay, when one of the two 
tau out of a Z-decay is identified. Therefore, taus are mostly left-handed, with some right- 
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Figure 2: i? distribution (defined as R = E n -/E Tj ) for NUHM scenarios and corresponding 
mSUGRA points. A cut R > 0.2 has been applied in each of these distribution. 
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handed admixtures. This shows up in the R-distribution for the corresponding mSUGRA 
points with a slight departure from that of the NUHM points. 

Designating the total number of events for 0.2 < R < 0.8 and R > 0.8 by 0\ and 2 
respectively, we find that the ratio r = Oi/0 2 is a rather effective discriminator between 
NUHM and a corresponding mSUGRA scenario yielding a similar number of same-sign ditau 
events. The values of this ratio for all the cases are listed in Table 3. For the NUHM points, 
this ratio turn out to be consistently larger than the corresponding mSUGRA points, which 
is expected from the R-distribution given in Figure 2. 

Lepton charge asymmetry: 

The results of charge asymmetry in the jet-lepton invariant mass distribution have been 
shown in Figure 3 and 4. For both the NUHM-1 and NUHM-3 benchmark points, gluinos and 
left-chiral squarks have closely spaced masses. Therefore, the hard jets are produced either 
in the decay qi —> qx\ 01 m Q~R ~^ QXii but not in whichever is allowed between g — > qqi,R 
or qi,R — > gq- This is due to small mass splitting between them; even if the gluinos are 
lighter than the left-chiral squarks, the decay chain q^ — > qx\ — > Qlt^ T ~^ Q^t^Xi is still 
the dominant source of the opposite sign same flavor dilepton signal, as the decay branching 
ratio of qi — > qg is very small (~ 2% or less) due to phase-space suppression. The branching 
fraction for q^ — > q\\ is ~ 32% and subsequently x\ decays into a l^l T pair with a decay 
branching fraction ranging from 21%-29%, while the sleptons decay into a lepton and the 
lightest neutralino with 100% branching ratio. 

It is clear from Figure 3a and 3b that both for NUHM-1 and NUHM-3 we get the 
desired charge asymmetry (which is negative for increasing rriji, since the lighter sleptons are 
dominantly left-chiral, and the leptons produced in decay are mostly back-to-back with 
the quark while the antileptons are in the same direction to that of the quark, hence re- 
distribution has larger population than m,ji+ distribution near the end-point of rriji invariant 
mass distribution. 

The situation is somewhat more complicated for the corresponding mSUGRA points. 
For mSUGRA- 1, the sleptons are heavier than the second lightest neutralino and the decay 
X 2 —> is suppressed. Here the main source of the flavor subtracted opposite sign same 
flavor dilepton signal is the two step q L — > qx\ — > q^l^Xi decay chain rather than the three 
step decay chain considered earlier. In this case X2 decays into a / ± / T Xi P arr y i a an off- 
shell slepton or Z. The sleptons are lighter than the second lightest neutralino and mostly 
dominated by the right-chiral component in mSUGRA-3. X2 follows its usual three step 
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Figure 4: a) rriji and b) A, vs rriji distribution for mSUGRA points. The mSUGRA-1 (top) and 
mSUGRA-3 (bottom) corresponds to NUHM-l,and NUHM-3 benchmark points, respectively. The 
event rates are predicted at an integrated luminosity of 10 fb" 1 . 
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decay chain. The expected positive charge asymmetry is visible for both the mSUGRA-1 
and mSUGRA-3 BP's. 

6 Summary and conclusion 

We have attempted a differentiation between mSUGRA and a scenario with non-universal 
Higgs masses. The extreme situation of large negative S, for which the characteristic features 
of the NUHM spectrum are most prominent, has been selected for this purpose, including 
three possible hierarchies among the masses of the lightest neutralino, the lighter stau and the 
tau-sneutrino. The primary channel of investigation being tau-rich, regions in the parameter 
spaces of both the scenarios, giving rise to similar ditau event rates, have been pitted against 
each other. 

In the same-sign ditau channel, we find that the ratio defined as R, the fraction of the 
energy carried by the charged pion in a jet produced in one-prong tau-decays, is a rather 
useful differentiator. Because of the dependence of R on the polarisation of the tau, one ends 
up having different numbers of events for the two cases in the regions R < 0.8 and R > 0.8. 
The ratios of these two event numbers, in turn, display a concentration in different regions, 
depending on whether it is NUHM or mSUGRA. 

We have further suggested the utilisation of signals involving leptons of the first two 
families, which are largely left chiral in NUHM. A bin-by-bin analysis of the of lepton- 
jet invariant masses exhibits a difference between the cases with negatively and positively 
charged leptons, whose general nature helps one distinguishing an NUHM scenario. 

If SUSY is indeed discovered at the LHC, one will certainly wish to run the machine 
with large integrated luminosity, so as to reveal the nature of the underlying scenario. One 
important question to ask in this context will be whether Higgs mass(es) have different high- 
scale origins compared to masses of the remaining scalars, namely, squarks and sleptons. A 
study in the line suggested here, based on the polarization study of tau as well as the first 
two family leptons, can be helpful in finding an answer to such a question. 

Acknowledgment : 

SB would like to thank Mihoko M. Nojiri and Theory Group, KEK and the Institute of 
Physics and Mathematics of the Universe for their hospitality while part of this work was 
being carried out. We also thank Nabanita Bhattacharya for her help during the preparation 
of this manuscript, and Atri Bhattacharya for computational assistance. This work was 



23 



partially supported by funding available from the Department of Atomic Energy, Government 
of India for the Regional Centre for Accelerator- based Particle Physics, Harish-Chandra 
Research Institute. Computational work for this study was partially carried out at the 
cluster computing facilities of KEK, Theory center and Harish-Chandra Research Institute 
(http : // cluster . mr i . ernet . in) . 

References 

[1] For reviews see for example, H. E. Haber and G. L. Kane, Phys. Rept. 117, 75 (1985). 

[2] S. Dawson, E. Eichten and C. Quigg, Phys. Rev. D 31, 1581 (1985); X. Tata, arXiv:hep- 
ph/9706307; M. E. Peskin, arXiv:0801.1928 [hep-ph]. 

[3] S. P. Martin, arXiv:hep-ph/9709356, and references therein. 

[4] A. Djouadi et al. [MSSM Working Group], arXiv:hep-ph/9901246. 

[5] A. H. Chamseddine, R. Arnowitt and P. Nath, Phys. Rev. Lett. 49, 970 (1982); 
R. Barbieri, S. Ferrara and C. A. Savoy, Phys. Lett. B 119, 343 (1982); L. J. Hall, 
J. Lykken and S. Weinberg, Phys. Rev. D 27, 2359 (1983); P. Nath, R. Arnowitt and 
A. H. Chamseddine, Nucl. Phys. B 227, 121 (1983); N. Ohta, Prog. Theor. Phys. 70, 
542 (1983). 

[6] M. Dine, W. Fischler, Nucl. Phys. B204, 346 (1982); J. R. Ellis, L. E. Ibanez, 
G. G. Ross, Nucl. Phys. B221, 29-67 (1983); C. Kounnas, A. B. Lahanas, 
D. V. Nanopoulos, M. Quiros, Nucl. Phys. B236, 438 (1984); A. A. Anselm and 
A. A. Johansen, "Susy GUT With Automatic Doublet - Triplet Hierarchy," Phys. 
Lett. B 200 (1988) 331; N. Polonsky and A. Pomarol, "GUT effects in the soft super- 
symmetry breaking terms," Phys. Rev. Lett. 73 (1994) 2292 [arXiv:hep-ph/9406224]; 
R. Hempfling, "Neutrino Masses and Mixing Angles in SUSY-GUT Theories with ex- 
plicit R-Parity Breaking," Nucl. Phys. B 478 (1996) 3 [arXiv:hep-ph/95 11288]. 

[7] J. R. Ellis, C. Kounnas and D. V. Nanopoulos, Nucl. Phys. B 247, 373 (1984); J. R. El- 
lis, K. Enqvist, D. V. Nanopoulos and K. Tamvakis, Phys. Lett. B 155, 381 (1985); 
M. Drees, Phys. Lett. B 158, 409 (1985); A. Corsetti and P. Nath, Phys. Rev. D 64, 
125010 (2001) [arXiv:hep-ph/0003186]; A. Corsetti and P. Nath, Phys. Rev. D 64, 



24 



125010 (2001) [arXiv:hep-ph/0003186]; S. Bhattacharya, A. Datta and B. Mukhopad- 
hyaya, JHEP 0710, 080 (2007) [arXiv:0708.2427 [hep-ph]]; S. Bhattacharya, A. Datta 
and B. Mukhopadhyaya, Phys. Rev. D 78, 115018 (2008) [arXiv:0809.2012 [hep-ph]]. 

[8] N. Chamoun, C. -S. Huang, C. Liu, X. -H. Wu, Nucl. Phys. B624, 81-94 (2002) [hep- 
ph/0110332]; K. Huitu, J. Laamanen, Phys. Rev. D79, 085009 (2009) [arXiv:0901.0668 
[hep-ph]]; S. P. Martin, Phys. Rev. D 79, 095019 (2009) [arXiv:0903.3568 [hep-ph]]; 
S. Bhattacharya, J. Chakrabortty, Phys. Rev. D81, 015007 (2010) [arXiv:0903.4196 
[hep-ph]]. 

[9] M. Drees, Phys. Lett. B 181 (1986) 279; J. S. Hagelin and S. Kelley, Nucl. Phys. B 342 
(1990) 95; Y. Kawamura, H. Murayama and M. Yamaguchi, Phys. Lett. B 324 (1994) 
52 [arXiv:hep-ph/9402254]; Y. Kawamura, H. Murayama and M. Yamaguchi, Phys. 
Rev. D 51 (1995) 1337 [arXiv:hep-ph/9406245]; A. Datta, A. Datta and M. K. Parida, 
Phys. Lett. B 431, 347 (1998) [arXiv:hep-ph/9801242]; A. Datta, A. Datta, M. Drees, 
D. P. Roy, Phys. Rev. D61, 055003 (2000), [arXiv:hep-ph/9907444]. 

[10] F. Gabbiani, E. Gabrielli, A. Masiero and L. Silvestrini, Nucl. Phys. B 477, 321 (1996) 
[arXiv:hep-ph/9604387]; M. Misiak, S. Pokorski and J. Rosiek, Adv. Ser. Direct. High 
Energy Phys. 15 (1998) 795 [arXiv:hep-ph/9703442]; J. Guasch and J. Sola, Nucl. 
Phys. B 562 (1999) 3 [arXiv:hep-ph/9906268]. 

[11] J. R. Ellis, K. A. Olive and Y. Santoso, Phys. Lett. B 539 (2002) 107, [arXiv:hep- 
ph/0204192]; J. R. Ellis, T. Falk, K. A. Olive and Y. Santoso, Nucl. Phys. B 652 
(2003) 259, [arXiv:hep-ph/0210205]. 

[12] H. Baer, A. Mustafayev, S. Profumo, A. Belyaev and X. Tata, JHEP 0507, 065 (2005), 
[arXiv:hep-ph/0504001]; J. R. Ellis, K. A. Olive and Y. Santoso, JHEP 0810 (2008) 
005, [arXiv:0807.3736 [hep-ph]]. 

[13] A. Katz and B. Tweedie, Phys. Rev. D 81 (2010) 035012, [arXiv:091 1.4132 [hep-ph]]; 
T. Figy, K. Rolbiecki and Y. Santoso, Phys. Rev. D 82, 075016 (2010), [arXiv: 1005.5136 
[hep-ph]]; 

[14] U. Chattopadhyay and D. Das, Phys. Rev. D 79, 035007 (2009) [arXiv:0809.4065 [hep- 
ph]]; S. Bhattacharya, U. Chattopadhyay, D. Choudhury, D. Das and B. Mukhopad- 
hyaya, Phys. Rev. D 81, 075009 (2010), [arXiv:0907.3428 [hep-ph]]. 

[15] M. M. Nojiri, Phys. Rev. D 51, 6281 (1995), [arXiv:hep-ph/9412374]. 



25 



[16] R. M. Godbole, M. Guchait and D. P. Roy, Phys. Lett. B 618, 193 (2005), [arXiv:hep- 
ph/0411306]; R. M. Godbole, M. Guchait and D. P. Roy, Phys. Rev. D 79, 095015 
(2009), [arXiv:0807.2390 [hep-ph]]; M. Guchait and D. P. Roy, arXiv:0808.0438 [hep- 
ph]. 

[17] S. Y. Choi, K. Hagiwara, Y. G. Kim, K. Mawatari and P. M. Zerwas, Phys. Lett. B 
648, 207 (2007), [arXiv:hep-ph/0612237]. 

[18] P. Richardson, JHEP 0111, 029 (2001), [hep-ph/0110108]. 

[19] A. J. Barr, Phys. Lett. B 596, 205 (2004), [arXiv:hep-ph/0405052]. 

[20] T. Goto, K. Kawagoe, M. M. Nojiri, Phys. Rev. D70, 075016 (2004), [hep-ph/0406317]; 
S. K. Mandal, M. Nojiri, M. Sudano and T. T. Yanagida, JHEP 1101, 131 (2011), 
[arXiv:1004.4164 [hep-ph]]. 

[21] A. Djouadi, J. L. Kneur and G. Moultaka, Comput. Phys. Commun. 176, 426 (2007), 
[arXiv:hep-ph/0211331]. 

[22] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667, 1 (2008). 

[23] A. Djouadi, M. Drees and J. L. Kneur, JHEP 0603, 033 (2006), [arXiv:hep- 
ph/0602001]. 

[24] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 192, 18 (2011), 
[arXiv:1001.4538 [astro-ph.CO]]. 

[25] J. R. Ellis, K. A. Olive, Y. Santoso and V. C. Spanos, Phys. Lett. B 588, 7 (2004), 
[arXiv:hep-ph/0312262]; J. L. Feng, S. Su and F. Takayama, Phys. Rev. D 70, 075019 
(2004), [arXiv:hep-ph/0404231]. 

[26] K. Hagiwara, A. D. Martin and D. Zeppenfeld, Phys. Lett. B 235, 198 (1990); 
B. K. Bullock, K. Hagiwara and A. D. Martin, Phys. Rev. Lett. 67, 3055 (1991). 

[27] S. Raychaudhuri and D. P. Roy, Phys. Rev. D 52, 1556 (1995), [arXiv:hep-ph/9503251]. 

[28] H. L. Lai et al. [CTEQ Collaboration], Eur. Phys. J. C 12, 375 (2000), [arXiv:hep- 
ph/9903282]. 



26 



[29] S. Jadach, J. H. Kuhn, Z. Was, Comput. Phys. Commun. 64, 275-299 (1990); Z. Was, 
P. Golonka, Nucl. Phys. Proc. Suppl. 144, 88-94 (2005), [hep-ph/0411377]; P. Golonka, 
B. Kersevan, T. Pierzchala, E. Richter-Was, Z. Was, M. Worek, Comput. Phys. Com- 
mun. 174, 818-835 (2006), [hep-ph/03 12240]. 

[30] T. Sjostrand, S. Mrenna and P. Skands, JHEP 0605, 026 (2006), [arXiv:hep- 
ph/0603175]. 

[31] M. Muhlleitner, A. Djouadi and Y. Mambrini, Comput. Phys. Commun. 168, 46 (2005) 
[arXiv:hep-ph/0311167]. 

[32] S. Biswas and B. Mukhopadhyaya, Phys. Rev. D 79, 115009 (2009), [arXiv:0902.4349 
[hep-ph]]. 

[33] G. Marchesini, B. R. Webber, G. Abbiendi, I. G. Knowles, M. H. Seymour and 
L. Stanco, Comput. Phys. Commun. 67, 465 (1992); G. Corcella et al, JHEP 0101, 
010 (2001), [arXiv:hep-ph/0011363]; S. Moretti, K. Odagiri, P. Richardson, M. H. Sey- 
mour and B. R. Webber, JHEP 0204, 028 (2002), [arXiv:hep-ph/0204123]. 

[34] F. E. Paige, S. D. Protopopescu, H. Baer and X. Tata, arXiv:hep-ph/0312045. 

[35] E. Richter-Was, arXiv:hep-ph/0207355. 

[36] The CMS Collaboration, CMS-TDR-8.1, CERN/LHCC 2006-001. 

[37] G. Aad et al. [The ATLAS Collaboration], [arXiv:0901.0512]. 



27 



